hep-ph/0209177 
September 2002 



Analytical Evaluation of Double BoxesQ 



V.A. SmirnovQ 



Nuclear Physics Institute of Moscow State University 
Moscow 119899, Russia 



Abstract 

Recent results on the analytical evaluation of double-box Feynman inte- 
grals and the corresponding methods of evaluation are briefly reviewed. 



1 Talk presented at the 36th Annual Winter School on Nuclear and Particle Physics (Repino, 
Russia, 25 February - 03 March 2002). 

2 E-mail: smirnov@theory.sinp.msu.ru. Supported by the Russian Foundation for Basic Research 
through project 01-02-16171, INTAS through grant 00-00313, and the Volkswagen Foundation, con- 
tract No. 1/77788. 



Introduction. 

Feynman diagrams with four external lines contribute to many important physical 
quantities. They are rather complicated mathematical objects because they are func- 
tions of multiple variables: internal masses, two independent Mandelstam variables 
s = (P1+P2) 2 , t = (P1+P3) 2 and squares of external momenta. The most complicated 
diagrams are the planar double box and non-planar (crossed) double box shown in 
Fig. 0. 




Figure 1: The planar and non-planar double boxes 

Almost all available analytical results correspond to the massless diagrams. Iron- 
ically, the first result for the massless double boxes was obtained in the most compli- 
cated case, where all the fours external legs are off-shell, i.e. p 2 7^ for all i — 1,2, 3, 4. 
This is the elegant analytical result for the scalar planar master double box, i.e. for 
all powers of propagators equal to one, obtained in M: 
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where 



C(xx,x 2 ,x 3 ) = - (6 [Li 4 (-pz) + Li 4 (-py)] 



y y 

+3 In - [Li 3 (-px) - Li 3 (-py)] + - In 2 - [Li 2 (-px) + U 2 (-py)] 
+~ ln 2 (px) \n 2 (py) + y ln(px) ln(py) + ^ In 2 | + ^- j , (2) 



AO, y) = J (I -x-y) 2 -4xy , p(x, y) = — , (3) 

v 1 — x — y + A(x, y) 

x = X1/X3, y = X2/X3, and Li n (^) is the polylogarithm 0. 

This diagram is convergent both in the ultraviolet and infrared sense. (For ex- 
ample, if one puts a dot on some line, an infrared divergence appears so that a 
regularization is needed.) To obtain this result, the authors have exploited the tech- 
nique of Feynman parameters and Mellin-Barnes (MB) representation. However, no 
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other results for the pure off-shell double boxes have been derived up to now so that 
the above result stays unique in the pure off-shell category. 

For massless double-box diagrams with at least one leg on the mass shell, i.e. pf = 
0, infrared and collinear divergences appear, so that one introduces a regularization 
which is usually chosen to be dimensional [H], with the space-time dimension d as a 
regularization parameter. One hardly believes that a regularized double-box diagram 
can be analytically evaluated for the general value of the regularization parameter 
e = (4 — d)/2, and the evaluation is usually performed in a Laurent expansion in e, 
typically, up to a finite part. 

The problem of the evaluation of Feynman integrals associated with a given graph 
according to some Feynman rules is usually decomposed into two parts: reduction of 
general Feynman integrals of this class to so-called master integrals (which cannot 
be simplified further) and the evaluation of these master integrals. A standard tool 
to solve the first part of this problem is the method of integration by parts (IBP) 
[|J when one writes down identities obtained by putting to zero various integrals of 
derivatives of the general integrand connected with the given graph and tries to solve 
a resulting system of equations to obtain recurrence relations that express Feynman 
integrals with general integer powers of the propagators through the master integrals. 

The purpose of this brief review is to characterize the status of the analytical 
evaluation of double-box diagrams and describe the corresponding techniques. We 
shall first deal with the pure on-shell case, i.e. pf = 0, i — 1,2,3,4, where this 
problem was completely solved during last three years, i.e. all the master integrals 
were calculated in expansion in e and a reduction procedure was developed both in 
the planar and non-planar case || ||, |7|, ||, |9|, [10], |Tl| . 

We also describe a method based on alpha (or Feynman) parameters and Mellin- 
Barnes (MB) representation used to calculate the master integrals, starting with a 
much simpler example of the on-shell box diagram. It turns out that the calculation 
of the basic master double-box diagram p[ based on a fivefold MB representation was 



far from being optimal. As it has been shown in |TT| it is possible to go through a 
fourfold MB representation. We shall describe here how this can be done. 

Then we turn to an intermediate massless case where one of the external legs is 
on-shell and the other three legs are on-shell, i.e. p\ = q 2 ^ 0, pf = 0, i — 2, 3, 4 
and list results of last two years where the problem of the evaluation was solved in 



this situation 12, 131. We then consider the massive on-shell case and describe the 



only available and recently obtained result for the master planar double box |TJ 



We conclude with a discussion of perspectives of the analytical evaluation of the 
double-box diagrams. 

The pure on-shell case: a box. 

To illustrate the method of MB representation let us evaluate the massless scalar 
on-shell box diagram of Fig. Q The Feynman integral can be written as 
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Figure 2: Box diagram. 
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F(s,t;d) j ( k2 + 2pi .k)(k 2 -2p 2 -k)k 2 (k + Pl +p 2 ) 2 ' (4) 

where usual prescriptions k 2 = k 2 + iO, etc. are implied. Using alpha or Feynman 
parameters we arrive at the following three-dimensional parametric integral: 

F(s, t- d) = m d ' 2 T(2 + e) f 1 f 1 f 1 d^i d6 d V V (l - r,) 

Jo Jo Jo 

x [- s £ l (i-Z l )(i- T] ) 2 -t&(l-&)ri 2 -i0]~ 2 ~ e . (5) 
Let us now use the MB representation 

1 1 1 r +ioa Y z 

dz-—r(X + z)T(-z), (6) 



(X + Y) x r(A) 2m J-ioo X x + Z 

where the contour of integration is chosen in the standard way: the poles with a 
r(. . . + z) dependence are to the left of the contour and the poles with a T(. . . — z) 
dependence are to the right of it. Representation (|^) is applied to the factor with 
square brackets in (^). As a result, the two terms in the square brackets in (^) raised 
to the power — 2 — e are replaced by a product of some powers of these terms (at the 
cost of introducing an extra integration), and, after evaluating parametric integrals 
in terms of gamma functions, we obtain 

m d l 2 1 r+ioo /f 

F(s,t:d) = , n , —, r / dz I - 

v ' ' ; (-s) 2+e r(-2e) 2ixi J-ioo \s 

xr(2 + e + z)T{l + z) 2 T(-l - e - z) 2 T{-z) , (7) 

where the same prescription for dealing with poles is implied, i.e. the poles of T(2 + 
e + z)Y{l + z) 2 are to the left of the integration contour and the poles of T(— 1 — e — 
z) 2 Y{— z) are to the right of it. In the case Ree < 0, we can choose this contour to 
be a straight line parallel to the imaginary axis, while in the case Ree > 0, a more 
complicated contour has to be chosen. Anyway, if we put e = 0, the integral becomes 
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ill-defined because the first pole with a T(. . . + z) dependence, i.e. at z = — 1, and the 
first pole with a T(. . . — z) dependence, i.e. at z = — 1 — e, glue together and there is 
no space to satisfy the prescriptions for the contours. 

This observation shows how the poles of the Feynman integral in e are generated, 
from the point of view of MB integrals. The next step is to pick up the pole in e 
by taking a residue, e.g. at the pole z = — 1 — e (with the minus sign, of course) 
and shifting the integration contour across this point. A resulting integral can be 
taken over the line at —1 <Re;z < 0. There is no gluing in this integral so that the 
integral can be safely expanded in a Taylor series in e. Every term of this expansion 
can be integrated by closing the integration contour to the right, taking residues at 
the points z = 0, 1, 2, . . ., and summing up a resulting series. Taking into account the 
terms up to e 1 and combining them with the value of the above residue we arrive at 
the following result: 



F(s,t\d) 



st l e 
Li 3 (-£)-21r^Li 2 (-t 



2 4tt 2 
- [ln(-s) + ln(-t)] - + 2 ln(-s) ln(-t) - — 



t 
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} + 0(e 2 ), (8) 



where 7e is the Euler constant. 



The pure on-shell case: the basic master diagram. 

Let us now consider the general on-shell planar double box diagram of Fig. 1, 
i.e. with a general irreducible numerator and powers of propagators. We choose 
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Figure 3: Planar double box diagram. 



this irreducible numerator and the routing of the external momenta as in |TTfl - see 
Fig. 3. For convenience, we consider the factor with (k+pi +P2+P3) 2 corresponding 
to the irreducible numerator as an extra propagator but, really, we are interested only 
in the non-positive integer values of a§. This general double box Feynman integral 
takes the form 



K(a!, 



j;s,t;e) 



d d kd d l 



(k 2 )^ [(k + pi) 2 )} a2 [(k + Pi+ p 2 ) 2 ] a3 
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[(k+ Pl + P2 +P3) 2 }- aS (q) 

[(I +P1+ P 2 ) 2 )] a4 [(I +P1+P2+ P3) 2 ] as {l 2 Y 6 [{k ~ l) 2 } a " ' 



where k and I are respectively loop momenta of the left and the right box. 

To resolve the singularity structure of Feynman integrals in e it is very useful 
to apply the MB representation (Q) that makes it possible to replace sums of terms 
raised to some power by their products in some powers, at the cost of introducing extra 
integrations. In []5j, |7|, [12| MB integrations were introduced directly in alpha/Feynman 
parametric integrals. It turns out more convenient to follow (as in JTTj]) the strategy 
of [0] and introduce, in a suitable way, MB integrations, first, after integration over 
one of the loop momenta, I, and complete this procedure after integration over the 
second loop momentum, k. After appropriate changes of variables we arrive at the 
following fourfold MB representation of @ (see also |TT||): 



K(a 1 ,...,a s ;s,t;e) -■ 



1^=2,4,5,6,7 r («i) r ( 4 - «4567 ~ 2t){~s)-^ 
I r~\~ioo (t\ W T{a 2 + w)T{- W )T{z 2 + zAT{z 3 + zA 



x / dwdz 2 dz 3 dz 4 

(2mr J-ioo 



x 



X 



X 



(27ri) 4 J-ioo \sj r(ai + z 3 + z 4 )T(a 3 + z 2 + z 4 ) 

T(4 - a 13 - 2a 28 - 2e + z 2 + z 3 )r(ai 238 - 2 + e + z 4 )r(a 7 + w - z A ) 

T(4 — a 46 — 2a 57 — 2e — 2to — z 2 — z 3 ) 
r(a 4 567 - 2 + e + to - z 4 )T(a s - z 2 - z 3 - z 4 )T{-w - z 2 - z 3 - z 4 ) 

r(4 - 01238 - 2e + w - z 4 )T(a 8 - w - ^ 2 - z 3 - z 4 ) 
T(a 5 + to + z 2 + z 3 + z 4 )T(2 - a 567 - e — w — z 2 ) 
r(4 - aw - 2a 28 - 2e + ^ 2 + z 3 ) 
xT(2 - a 457 - e - to - z 3 )T(2 - a 128 - e + £ 2 )r(2 - a 238 - e + z 3 ) 
xT(4 - a 46 - 2a 57 - 2e - 2io - z 2 - z 3 ) , (10) 

where a 4 56 7 = a 4 + a 5 + a6 + a 7 , a 13 = a 4 + a 3 , etc., and integration contours are chosen 
in the standard way. 

In the case of the master double box, we set = 1 for i — 1, 2, . . . , 7 and a 8 = 
and obtain 

^°)( S ,t; e )=^(l,...,l,0; S ,t;e) 

( i7rd/2 ) 1 /■+*», , , , ft\ w T(l + w)T(-w) 

1 dwdz 2 dz 3 dz 4 



x 



r(-2e)(-s) 3 + 2e (2vri) 4 7-ioo W r(l - 2e + to - z 4 ) 

r(2 + e + to - z 4 )r(-l - e - to - z 2 )r(-l - e - to - ^ 3 ) 



r(i + 2 2 + z 4 )r(i + 2 3 + z 4 ) 

xT(l + to + z 2 + z 3 + z 4 )r(l + e + 2 4 )r(z 2 + z 4 )r(2; 3 + z A ) 

xr(-e + zjj)r(-e + z 3 )r(l + to - z 4 )T(-z 2 -z 3 - z 4 ) . (11) 

Observe that, because of the presence of the factor r(— 2e) in the denominator, we 
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are forced to take some residue in order to arrive at a non-zero result at e = 0, so 
that the integral is effectively threefold. 

Of course, the resolution of singularities in e in such a multi-dimensional MB 
integral is much more complicated than in the one- dimensional case. The poles in e 
are not visible at once, at a first integration over one of the MB variables. However, 
the rule for finding a mechanism of the generation of poles is just a straightforward 
generalization of the rule used in the previous one-loop example, where we saw that 
the product of T(— 1 — e — z)T(l + z) generated the pole of the type T(— e) (this is 
nothing but the value of one of these gamma function at the pole of the other gamma 
function). Suppose we start from the integration over one of the variables, z. We 
analyze various products T(a + z)T(b — z), where a and b depend on the rest of the 
variables, with the understanding that this integration generates a pole of the type 
T(a + b). This means that any contour of one the next integrations should be chosen 
according to this dependence. We continue this analysis with various integrations at 
the second step, etc. 

Here is an example of the procedure of generating poles in the integral ([IT]). The 
product r(— 1 — e — w — z 2 )Y{— e + z 2 ) generates, due to the integration over z 2 , a pole 
of the type T(— 1 — 2e — w). Then the product of this gamma function with T(l + w) 
generates a pole of the type T(2e) due to the integration in w. 

After such preliminary analysis we conclude that the key gamma functions that are 
responsible for the generation of poles in e are r(— e + z 2 ), T(—e+z 2 ) and r(l+io — z 4 ). 
This gives a hint for the construction of a complete procedure of the resolution of 
the singularities in e, with the goal to decompose the given integral into pieces where 
the Laurent expansion of the integrand in e becomes possible. One can proceed as 
follows. 

We first take care of the gamma functions T(— e + z 2 ) and T(— e + 23), i.e. take 
residues at z 2 = e and z 3 = e and shift contours across these poles. As a result, flllD is 
decomposed as K = K\\ + K w + K m + K 00 , where K\\ corresponds to taking the two 
residues, K 00 is defined by the same expression ( |TT|) but with both first poles of the 
selected two gamma functions treated in the opposite way, and the two intermediate 
contributions are defined by taking one of the residues and changing the nature of 
the first pole of the other gamma function. 

The contribution K\\ takes the form 



The contributions K w and K m are equal to each other because of the symmetrical 
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xF(-w)F(l +w- z 4 )T(2 + 6 + w- z 4 )T(e + z 4 ) 2 T(-2e - z A ) 
r(l + 2e + w + z 4 ) 



(12) 



r(l - 2e + w - z 4 )r(l + e + z 4 ) ' 
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dependence of the integrand on z 2 and z 3 . We have 



Km = -T(-2e)(-s)^(2nirJ- ioo dwdz ^[- s ) r(l + wn-l-2e-w) 
xr(-tu)r(-l -e-w- z 2 )T(-e + z 2 )T{l +w- z 4 )T(2 + e + w - z 4 ) 
T(e + z 4 )T(z 2 + z 4 )Y(-e — z% — z 4 )T(l + e + w + z 2 + z 4 ) 



T(l -2e + w-z 4 )T(l + z 2 + z 4 ) 



(13) 



where the first pole of T(— e + z 2 ) is of the opposite nature. 

For all these contributions, further decompositions are necessary. One can proceed 
as follows. 

{Ku) Take minus residue at w = — 1 — 2e. The resulting one- dimensional integral 
is calculated by taking care of the first poles of T(z 4 ) and T(z 4 + e). One obtains either 
an explicit expression in gamma and ip functions and a one- dimensional MB integral. 
In the integral, where the first pole of T(— 1 — 2e — w) is of the opposite nature, 
take care of the first pole of T(z 4 + e). One obtains a one-dimensional MB integral 
over w which is calculated by expanding in e. The integral, where the first poles of 
T(— 1 — 2e — w) and T(z 4 + e) are of the opposite nature, does not contribute because 
one can expand it in e and obtain a zero result due to T(— 2e) in the denominator. 

{Kqi) Take minus residue at w = — 1 — 2e and consecutively take care of the first 
poles of the gamma functions T(z 4 + e), T(z 2 + z 4 ) and T(z 2 + z 4 — e) in the resulting 
integral. One obtains one- dimensional MB integrals and a two-dimensional integral 
in z 2 and z 4 which is calculated by use of the second Barnes lemma (see below). In 
the integral with the pole of T(— 1 — 2e — w) of the opposite nature, one takes care 
of the first pole of T(z 2 + z 4 ), i.e. takes a residue at z 4 = —z 2 (and then takes care 
of the first pole of T(— 1 — e — w — z 2 )) and considers an integral with the opposite 
nature of this pole (with a zero contribution). 

(Kqq) Take care of the first poles of T(— 1 — e — w — z 2 ) and T(— 1 — e — w — z 3 ). 
The only non-zero contribution arises when taking both residues. 

As a result we obtain either explicit expression in terms of gamma functions 
and their derivatives, or one- dimensional integrals over straight lines parallel to the 
imaginary axis of ratios of gamma functions which can be of two types: integrals 
over w or some z-variable. The integrals over w can be calculated by closing contour 
to the right, taking residues at the points w = 0, 1, 2, . . . and summing up resulting 



series, with the help of the table of formulae presented in |fL5 |. The one-dimensional 
MB integrals in Zi can be calculated with the help of the first and the second Barnes 
lemmas 



— / dw T(Ai + w)T(X 2 + w)T(X 3 - w)T(X 4 - w) 

_ r_(Ai + A 3 )r(Ai + A 4 )r(A 2 + A 3 )r(A 2 + a 4 ) 

T(Ai + A 2 + A 3 + A 4 ) 



(14) 
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i /-n°° r(Ai + w)r(A 2 + w)r(A 3 + w)r(A 4 - w)r(A 5 - w) 

27TC 7-ioo T(Ai + A2 + A3 + A4 + A5 + w) 

_ r(Ai + A 4 )r(A 2 + A 4 )r(A 3 + A 4 )r(Ax + A 5 )r(A 2 + A 5 )r(A 3 + a 5 ) 

T(Ai + A 2 + A 4 + A 5 )r(Ai + A 3 + A 4 + A 5 )r(A 2 + A 3 + A 4 + A 5 ) 
and their corollaries. These are two typical examples of such corollaries: 



(15) 
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2m 



dw r(Ai + w)r(A 2 + w) 2 T(-\ 2 - w)T(X 3 - w) 
-ioo r(Ai + A 2 + A3 + w) 

r(Ai - A 2 )r(A 2 + A 3 ) W (Ai + A 3 ) - f (A 2 + A 3 )] 



(16) 



where the nature of the pole at w — - 
are treated in the standard way, and 



r(Ai + a 3 ) 

•A 2 is determined by r(A 2 + u>) while other poles 



1 ,-1/2+ioc 

— / dwT(l + w)T(w)T(-w)T(-l - wWl + wY 

2m 7-1/2-ioo 



ikl + 67eC(3) + ?L . 
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(17) 



Here ((z) is the Riemann zeta function. 

Collecting all the contributions, one reproduces the result of ||: 



where 
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+ - In 4 x + 6vr 2 In 2 x - —((3) lnx + — tt 4 
3 3 s v ; 30 

- 2 Li 3 (-x) - 2 In x Li 2 (-x) - (in 2 x + 7r 2 ) ln(l + x) 

-4 [^(-x) -lnx£ 1)2 (-x)] + 44Li 4 (-x) 
-4 [ln(l + x) + 6 lnx] Li 3 (-x) 

+2 (^ln 2 x + 2 In x ln(l + x) + y tt 2 ^) Li 2 (-x) 
+ (ln 2 x + 7r 2 ) ln 2 (l + x) 



41n 3 x + 5vr 2 lnx-6C(3) ln(l + x) + 0(e) 



(19) 
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where, in addition to the polylogarithms, one meets also the generalized polyloga- 



rithms 16 



S (z) - { ~ ir+b " f 1 ^'Wa-^h , (20) 
Sa > b{z) - (a-l)!6! Jo 1 dt ■ {20) 

This result is in agreement with the leading behaviour in the (Regge) limit t/s — > 
obtained by use of the strategy of expansion by regions JC], [18], [L{|. It turns out that 
the (h-h), (lc-lc) and (2c-2c) contributions are the only non-zero contributions. (See 
|18| , pTQfl for definitions of the hard (h) and collinear ((lc) and (2c)) regions and the 
corresponding contributions.) Keeping the two leading powers in x we have || 

4 51nx ( , 9 5 9 \ 1 

= - ? + ^-(2^ 2 X--7r 2 J- 

/2, , 11 2l 65 A/ A 1 

4 , ,* 9,0 88 , , . , 29 / 

+ -ln 4 x + 6vr 2 ln 2 x £(3) lnx+ — vr 4 

3 3 s v ; 30 

+2x (ln 2 x - 21nx + vr 2 + 2) 

-i {4 In 3 x + 3 In 2 x + (5vr 2 - 36) lnx + 2 [33 + 5tt 2 - 3C(3)]}) 
+0(x 2 ln 3 x,e) . (21) 



The pure on-shell case: the non-planar case and reduction to master 
diagrams. 

The basic master non-planar on-shell massless double-box diagram were calcu- 
lated in by the same method of Feynman parameters and MB representation, with 
the only qualification that the initial integration contours were chosen, at appropriate 
values of e, as straight lines parallel to imaginary axis. Then the procedure of taking 
residues and shifting contours was applied, with the requirement to keep this prop- 
erty.0 It turns out that it is natural to consider non-planar double boxes as functions 
of the three Mandelstam variables s,t and u = (pi + p^) 2 not necessarily restricted 
by the physical condition s + t + u = which does not simplify the result. 

Reduction procedures for the evaluation of general double-box diagrams, with 
arbitrary numerators and integer powers of the propagators were developed in 
in the planar case and in |J in the non-planar case. Both reduction procedures 
were based on solving recurrence relations following from IBP [|J. In the non-planar 
case, the so-called Lorentz invariance (LI) identities were used, in addition to the 
IBP relations. In 0, the first of the two most complicated master integrals in- 
volved is K^> = K(l, . . . , 1,0) considered above. As a second complicated master 

3 I think, this is a matter of taste what variant of resolution of the singularities in e to apply. 
Anyway, I have confirmed the result of using the strategy described above. 
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integral, the authors of B have chosen the diagram with a dot on the central line, 



i.e. K(l, 1, 1, 1, 1, 2, 1, 0). As was pointed out later pl |, in practical calculations one 
runs into a linear combination of these two master integrals with a coefficient 1/e, 
so that a problem has arisen because the calculation of the master integrals in one 
more order in e looked rather nasty. Two solutions of this problem have immediately 



appeared. In |10 |, the authors have calculated this very combination of the master 
integrals, while in JlTJ another choice of the master integrals has been made: instead 
of K(l, 1, 1, 1, 1, 2, 1, 0), the authors have taken the integral K(l, 1, 1, 1, 1, 1, 1, —1) as 
the second complicated master integral. This was a more successful choice because, 
according to the calculational experience, no negative powers of e occur as coefficients 
at these two master integrals. 

These analytical algorithms have been successfully applied to the evaluation of 



two-loop virtual corrections to various scattering processes |22| in the zero-mass ap- 
proximation. 

One leg off-shell. 

In the case, where one of the external legs is on-shell, p\ ^ 0, pf = 0, i = 2, 3, 4, 
the planar double box and one of two possible non-planar double-box diagram with all 



powers of propagators equal to one were analytically calculated in [[12], as functions 
of the Mandelstam variables s and t and the non-zero external momentum squared 
p\. Explicit results were expressed through (generalized) polylogarithms, up to the 
fourth order, dependent on rational combinations of p\ , s and t, and a one- and (in 
the non-planar case) two-dimensional integrals with simple integrands. To do this, 
the method based on MB integrals described above was applied. 

These and other master planar and non-planar double boxes with one leg off- 



shell were evaluated in [|13| with the help of the method of differential equations |2c 
The corresponding results are expressed through so-called two-dimensional harmonic 
polylogarithms which generalize harmonic polylogarithms introduced in ||24|| . 

A reduction procedure that provides the possibility to express any given Feynman 
integral to the master integrals was also developed in |T3|] . It is based on the observa- 



tion that, when increasing the total dimension of the denominator and numerator in 
Feynman integrals associated with the given graph, the total number of IBP and LI 
equations grows faster^ than the number of independent Feynman integrals (labeled 
by the powers of propagators and powers of independent scalar products in the numer- 
ators). Therefore this system of equations sooner or later becomes overconstrained, 
and one obtains the possibility to perform a reduction to master integrals. 

In fact, this strategy is quite general and can be applied to any Feynman graph. 
However its implementation for any concrete graph looks rather non-trivial. In partic- 
ular, solving an overconstrained system of several thousand equations is a challenging 
technical task. Nevertheless, this method was successfully applied |26| to the Feyn- 
man integrals with one leg off-shell contributing to the process e + e~ — > 3jets. 



4 This fact was first pointed out by Laporta and used in p5| 
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The massive on-shell master double box. 

We now turn to the massive on-shell case, i.e. pf = m 2 , i — 1,2, 3, 4, and consider 
the general on-shell planar double box diagram of Fig. |], i.e. with general irreducible 
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Figure 4: Planar double box diagram. Solid and dashed lines denote massive (with 
the mass m) and massless propagators, respectively. 

numerator and powers of propagators. The irreducible numerator, the numbering of 
lines and the routing of the external momenta are chosen as in the massless case of 
(|9|) (see Fig. [3]). This general double box Feynman integral takes the form 

d d kd d l 



B(a u 



s,t,m ; e) 



(F 



m 



2\ai 



[(k + Pi) 2 )} a2 [(k + Pl + p 2 ) 



m 



21a 3 



X- 



[(k + P1+P2+ Ps) 



-u 8 



[(/ + Pi + P2) 2 - m 2 )} ai [{I + P1+P2+ Pz) 2 ] a5 (I 2 - m 2 ) a6 [{k - l) 2 } a " ' (22) 

To arrive at a MB representation (with, possibly, minimal number of MB integra- 
tions) one straightforwardly generalizes the above procedure for the massless case by 
introducing two extra MB integrations, when separating terms with m 2 after each of 
the integrations over the loop momenta, and after appropriate changes of variables 



one obtains the following sixfold MB representation of (|22|) : 



B(a u ...,a 8 ;s,t,m 2 ;e 



(m d l 2 ) (-l) a 



x 



(2vn) 



dw Y[ d z j 

3=1 



^=2,4,5,6,7 r(a i )r(4 - a 45 67 - 2e) (s) a -^ 
'mT\ zl+Z5 ny T(a 2 + w)T(-w)T(z 2 + z 4 )T(z 3 + z A ) 



T(a 1 + z 3 + z 4 )T(a 3 + z 2 + z 4 ) 



x 



T(4 - a 13 - 2a 28 -2e + z 2 + Z3)r(ai238 - 2 + e + z 4 + z 5 )T(a 7 + w - z A ) 

T(4 — a 46 — 2a 57 — 2e — 2w — 2z\ — z 2 — z 3 ) 
r(a.4567 - 2 + e + w + zx- 2 4 )r(a 8 - z 2 — z 3 - z A )T(-w - z 2 - z 3 — z 4 ) 

T(4 - ai238 - 2e + w - -24)r(a 8 - w - z 2 - z 3 - z A ) 
T(a 5 + w + z 2 + z 3 + z 4 )r(2 - a 567 - e — w — z x — z 2 ) 
r(4 - a 13 - 2a 28 - 2e + z 2 + z 3 - 2z 5 ) 
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xT(2 - a 457 -e-w-zi- z 3 )T(2 - a 12S - e + z 2 - z 5 )T(2 - a 238 - e + z 3 - z 5 ) 
xr(4 - a 46 - 2a 57 - 2e - 2u> - z 2 - z 3 )T (- Zl )T (- z 5 ) . (23) 

In the case of the master double box, we set dj = 1 for i = 1, 2, . . . , 7 and a 8 = 
and obtain a massive generalization of (|IlD 

B (0) (s,t,m 2 ;e) = B(l, . . . , 1, 0; s,t,m 2 ; e) 



(m d / 2 )" 



A /m 2 V 1+ * 6 r(i + w)r(-w) 



r(-2e)(-s) 3+2e (2^)6^00 du, n d ^ (^_ s y r(l - 2e + w - z A ) 

T(2 + e + w + Z\ — Z/±)T{— 1 — e — to — Zi — z 2 )r(— 1 — e — u> — 24 — z^)Y{—Z\) 

r(l + z 2 + z 4 )T(l + z 3 + z 4 )T(-2e + z 2 + z 3 - 2z 5 ) 
T(-e + z 2 - z 5 )V(-e + z 3 - z 5 )T(l + £ + 2:4 + z 5 )V(-z 5 )V(-2e + z 2 + z 3 ) 



x 



T(-2 -2e-2w- 2z x - z 2 - z 3 ) 
xT(-2 - 2e - 2w - z 2 - z 3 )T(l + w + z 2 + z 3 + z 4 )T{z 2 + z i )T{z 3 + z 4 ) 
xT(l + w-z A )T{-z 2 -z 3 -z A ). (24) 



The resolution of singularities in e is performed |T4| also in the standard way (see 
H [7J 0, |T2[ and the examples above) and reduces to shifting contours and taking 
residues. The final result takes the following form: 



B^(s,t,m 2 ;e) 



i7T d/2 e - lE e\ x 2 r 



s 2 (-ty+ 



2c 



^ + ^ + b 01 (x) + b 02 (x, y) + 0{e) 



b i(x) 



where x = — Am 2 / 's, y = l/\Jl — Am? /t 

b 2 (x) 



2(m x -p x ) 2 , 

T . fl-x\ T . fl + x\ /-2x\ ( 2x 

Ll 3 ~ ^~ + Ll 3 — — + Ll 3 + Ll 3 



(25) 
(26) 



1 — x 



1 + x 



+A(m x - p.. 



Lie 



1 — x 



Lie 



-2x 



2 J VI - x 

2 , /o /qn„3 , ai r, 1 J2 



- (4/3)m^, + Am 2 xPl 



-6m x p 2 x + (2/3)p* + Al 2 (m xPx + p 2 x ) - 2l 2 (m x + 3p x ) 
-(tt 2 /3)(4/ 2 - m x - 3p x ) + (8/3)$ + 14Cs , 



(27) 



-tt{m x -p x 
, 2x 

-L13 



Li 3 (x) - Li 3 (-x) - Li 3 ( — ^) + Li 3 ( ' 



+ Li s 



-2x 



1 + VI — x 

8Li 2 (x) Li 2 (— x) + 16Li 2 



+ 4 
1 — x 



Li 2 (x) 2 + Li 2 {-xf + 4Li 2 1 ' ~ ' 



(Li 2 (x) - Li 2 (-x)) 



-(4/3) [tt 2 - 6/ 2 + 3m 2 + 6m,.(2/ 2 - 2l x - p x ) + I2l x p x - 3p 2 ](Li 2 (x) - Li 2 {-x)) 



(8/3) [7r 2 - 6/2 + 6/^ - 6m x (l x +p x - 2/ 2 )]Li 2 



1 — x 
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+8(m x - p x ) 



2/ 2 )Li 2 



2x 



+ 2(l x -m x + l 2 )U 2 



' -2x s 
.1 — x, 



-S(m x -p x )(2l x -p x - 5m x + Al 2 )(-m x p x + l 2 (m x + p x ) - l\ + vr 2 /6) 
-(20/3)m* + (164/3)m*p x - 40m 2 p 2 _ (4/3)m^ - (8/3)p x 
+8m x l x (m 2 x - 3m x p x + 2p 2 x ) 

-Al 2 (7m 3 x + 21m 2 x p x - Am x l x p x - 23m x p 2 x + Al x p 2 x - p 3 x ) 
-7r 2 ((17/3)m 2 - {A/3)mJ x - 2m xPx + (4/3)Z x p x - (7/3)p 2 ) 
+/2(84m 2 - 8)71,4 - 16m x p x + 8Z x p x - 44p 2 ) 

— (8/3)Z 2 (6Z| - 7r 2 )(3m x - 2 Pa; ) - (4/3)tt 2 / 2 + &\ + vr 4 /9 , (28) 



and 



b 2(x,y) = 2{p x -m x )\A 



Li a 



1 — x 



Lis 



1 + x 



+ Li, 



-Li a 



Lis 



+2 
-Li 3 



Li 3 



(1 +x)y 
1-xy / 

f (i+ x )(i- y y 

V 2(1 -xy) , 
(l-x)(l-y)' 



-(i -^)y x 

1 + xy 



Li a 



(l_-x)y N 
1 - xy t 

' (l + x)y 
1 + xy 



Li- 



2(1 + xy) 



+ Li 3 



/ (l-x)(l + y) ' 
''[ 2(1 -ay) , 
(l+x)(l+yV " 
2(1 + xy) ( 



+2(m y + p y - m X2/ - p X2/ ) 



2Li 2 (x) - 2Li 2 (-x) + Li 2 



-2x 
1 — x 



Li? 



+A(m xy - p xy )(U 2 (-y) - Li 2 (y)) - A(m x + p x - 2/ 2 )Li 2 ^ 



2x 

1 + x 
1 — x 



-A(m xy - p xy )U 2 



+4(p x + l y - m xy )Li 2 



i-y 



-(l + x)y N 



4(m x + /j, - m xy )Li 2 



+4(p x + Z w - p^ y )Li 2 



1 — xy 
(l + x)y 



- A(m x + l y -p xy )Li 2 



' (1 -x)y \ 
1-xy J 

-(1 - x)y 
1 + xy 



1 + xy 



+2(m x + p x + m y + p y - 2m xy - 2/ 2 )Li 2 
+2(m x +p x + m y +p y - 2p xy - 2/ 2 )Li 2 



f (i-x)(i+ y y 

V 2(1 -xy) 
(l-x)(l-y) \ 
2(1 + xy) y 



+2p 2 (m J/ +p y - m xy - p xy ) + 2p x (2(m y l y + m y p y + Z y p y ) 
+m X y(-m y - 2/ y - 3pj, + 3m xy ) + p xy (-3m y - 2l y - p y + 3p xy )) 



+2m x (2p x + m v - 2l y + p y )( 
+2p y (2m 2 + p 2 ) + m 2 {2p 



Py 



m 



xy 



Vxy) -P 2 y(m xy +P X y) 



'll '"j-!/ Vxy) 
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+2m y (p 2 y + m 2 xy + 2p 2 xy - p y (3m xy + p xy )) - 2(m 3 xy + p\ ; 
+2l 2 ((Am y + Ap y - 3m xy )m xy + (2m y + 2p y - 3p xy )p xy 



-2{p x + 2m x ) (m y + p y - m xy - p xy ) -m\- Am y p y - p 2 y ) 
+2l 2 2 (3{m y + p y ) - 2(2m xy + p xy )) 
-(ir 2 /3)(m y +p y - 8m xy + Qp xy )} ■ 



(29) 



Here the following abbreviations are used: £3 = C(3), l z = lnz for z = x,y,2, p z 
ln(l + z) and m z = ln(l — z) for z — x,y, xy. 



The result (^)-(|29|) is in agreement with the leading power behaviour in the (Su- 
dakov) limit of the fixed-angle scattering, m 2 -C \s\, \t\ 



This asymptotics is obtained 
by use of the strategy of expansion by regions [|17], 0, |l9| . The structure of regions 
is very rich. The following family of seventeen regions participates here: 

(h-h), (lc-h), . . . , (4c-h), (lc-lc), . . . , (4c-4c), 
(lc-3c), (2c-4c), (lc-4c), (2c-3c), 
(luc-2c), (2uc-lc), (3uc-4c), (4uc-3c). 

Here h denotes hard, c - collinear and uc - ultracollinear regions for the two loop 
momenta. (See [113] and Chapter 8 of [H5 for definitions of these regions.) In par- 



ticular, the (h-h) contribution is nothing but the massless on-shell double box fliTf). 



Evaluating and summing up all the contributions, we obtain |14 



B(°\s,t,m 2 -e) 



/ i7r d/2 e -7EeY 



l+2e 




(2/3)L 3 + (tt 2 /3)L + 2C 3 



+ 



4Li, 



-t/s) 



-(2/3)L 4 + 2 ln(t/s)L 3 - 2(ln 2 (t/s) + An 2 /3)L 2 
4 ln(t/s)Li 2 (-t/s) + (2/3) \n 3 (t/s) - 2 ln(l + t/s) \n 2 (t/s) 



+ (8tt 2 /3) \n(t/s) - 2tt 2 ln(l + t/s) + IOC3] L + tt 4 /36} + 0(m 2 L 3 , e) , (30) 

where L = ln(— m 2 /s). This asymptotic behaviour is reproduced when one starts 
from result 

It is interesting to note that, in the above result, there are no so-called two- 



dimensional harmonic polylogarithms |24| which have turned out to be adequate 



functions to express results for the double boxes with one leg off-shell fT3 |. 



Perspectives. 

It is believed that sooner or later we shall achieve the limit in the process of 
analytical evaluation of Feynman integrals so that we shall be forced to proceed only 



numerically. (See, e.g., the first paper of []27 where this point of view has been 
emphasized.) The progress in the field of numerical evaluation of Feynman integrals 
was rather visible last years. Several new powerful numerical methods have been 
developed — see, e.g., |29l [30|, p7fl . They have been applied in practice and, at the 



same time, served for crucial checks of analytical results. For example, the method 
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of |28| applicable to Feynman integrals with severe ultraviolet, infrared and collinear 
divergences was successfully applied to check a lot of results for the double-box master 
integrals presented and/or discussed above. 

However the dramatic progress in the field of analytical evaluation of Feynman 
integrals shows that we have not yet exhausted our abilities in this direction. Let 
us take the problem of the analytical evaluation of the on-shell double-box diagrams 
with p 2 = m 2 , i = l,2,3,4 as an example. The calculational experience, in particular 
obtained when evaluating four-point diagrams, tell us that if such master integrals 
can be evaluated, the problem can be also completely solved, after evaluating other 
master integrals and constructing a recursive procedure that expresses any given Feyn- 
man integral with general numerators and integer powers of propagators through the 
master integrals. Therefore the explicit analytical result (|25|)-(|29|) can be considered 
as a kind of existence theorem, in the sense that it strongly indicates the possibility 
to analytically compute a general Feynman integral in this class (and apply these 
results to various scattering processes in two loops without putting masses to zero). 
To calculate the master integrals one can apply the technique of MB integration de- 
scribed above. To construct appropriate recursive algorithms both in the planar and 
non-planar case one can use recently developed methods based on shifting dimension 
20| and differential equations |13[ as well a method based on non-recursive solutions 



of recurrence relations |31] (successfully applied in practice |E 

One can also hope that new analytical results can be obtained for many other 
classes of Feynman integrals depending on two and three scales. In particular, the 
analytical evaluation of any two-loop two-scale Feynman integral with two, three and 
four legs looks quite possible. 
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